GLAST sensitivity to Point Sources of Dark Matter Annihilation 
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We study the prospects for detecting gamma-rays from point sources of Dark Matter annihilation 
with the space satellite GLAST. We simulate the instrument response to the gamma-ray spectrum 
arising from the annihilation of common Dark Matter candidates, and derive full-sky sensitivity 
maps for the detection of point sources and for the identification of the Dark Matter (as opposed 
to astrophysical) origin of the gamma-ray emission. These maps represent a powerful tool to assess 
the detectability of point sources, i.e. sources with angular size smaller than the angular resolution 
of GLAST, AO ~ 0.1°, in any DM scenario. As an example, we apply the obtained results to 
the so-called mini-spikes scenario, where the annihilation signal originates from large Dark Matter 
overdensities around Intermediate Mass Black Holes. We find that if these objects exist in the 
Galaxy, not only GLAST should be able to detect them over a timescale as short as 2 months, but 
in many cases it should be possible to determine with good accuracy the mass of the annihilating 
Dark Matter particles, while null searches would place stringent constraints on this scenario. 

PACS numbers: 



I. INTRODUCTION 

In the cosmological concordance model that has 
emerged over the past few years, a considerable frac- 
tion of the total energy density today, Qcdm = 0.19, 
consists of non-baryonic, collisionless and dissipation-free 
(i.e. cold) matter [l|]. To unveil the - so far completely 
unknown - nature of this dark matter (DM) is one of the 
most outstanding challenges for cosmology and astropar- 
ticle physics today. The list of proposed DM candidates 
is long, ranging from modified theories of gravity, that 
would effectively mimic a large non-relativistic compo- 
nent in the total energy content of the universe, to a 
whole zoo of speculative new particles that the DM may 
consist of. While the former approach may, in fact, be 
extremely successful to describe certain isolated phenom- 
ena, like the flattening of galactic rotation curves , it is 
notoriously difficult to reconcile with the whole range of 
accessible observations; in the following, we will therefore 
restrict ourselves to the latter possibility (see Refs. 0, Q 
for recent reviews on particle DM). 

Search strategies for DM particles (going beyond 
merely testing their gravitational influence) can be 
grouped into two categories. In direct detection experi- 
ments, one tries to trace these particles by looking for the 
recoil energy they would transfer during scattering events 
with the atoms of the detector material; null searches 
start to place non-trivial bounds on the parameter space 
of the underlying models (see, e.g., Q and references 
therein). Alternatively, one can use indirect detection 
techniques, making use of the fact that DM particles will 
generally pair-annihilate in regions of enhanced DM den- 
sities; the decay products may then be revealed as exotic 
contributions to astrophysical fluxes of gamma-ray, neu- 
trino and anti-matter. 

The space satellite GLAST is expected to play a crucial 
role in indirect DM searches, thanks both to its ability to 



perform observations at energy scales comparable to the 
mass of common DM candidates and to its potential of 
making deep full-sky maps in gamma-rays, thanks to its 
large (~ 2.4 sr) field-of-view Q. In the first part of this 
paper, we study the prospects for detecting point sources 
of DM annihilation with GLAST, and assess the possi- 
bility to discriminate them from ordinary astrophysical 
sources. As a result of this analysis, we obtain full-sky 
sensitivity maps for the detection of sources above the 
diffuse background, as well as for the identification of 
DM annihilation sources, i.e. for the detection of sources 
that can be discriminated against ordinary astrophysical 
sources. 

This analysis can be applied to any astrophysical sce- 
nario that predicts the existence of point-like sources of 
gamma-rays from DM annihilations, where "point-like" 
in this context means sources with an angular size smaller 
than the GLAST angular resolution, A9 ~ 0.1° (see be- 
low for further details). In the second part of the pa- 
per, in order to show the effectiveness of this approach, 
we apply our results to a specific astrophysical scenario, 
where gamma-rays originate from DM mini-spikes. In 
fact, it has been recently pointed out that intermedi- 
ate mass black holes (IMBHs) are a particularly promis- 
ing place to look for DM self-annihilations, since, under 
certain assumptions about their formation mechanism, 
they will be surrounded by regions of highly enhanced 
DM density, called mini-spikes (34[. As a consequence, 
they are expected to appear as bright sources of gamma 
rays, thereby providing a new possible means of indirect 
DM detection. Using the sensitivity maps obtained in 
the first part of the present paper, we study in detail 
the prospects for detecting these objects, and show that 
GLAST can find a large fraction of the mini-spike pop- 
ulation in the Milky Way over a time scale as short as 
2 months, or set strong constraints on their existence in 
case of null searches. Finally we show that, in most cases, 
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one can not only expect a clear discrimination from ordi- 
nary astrophysical sources, but also obtain a satisfactory 
determination of the DM mass as well as hints on its 
nature, through the determination of the leading annihi- 
lation channel. If the mass of the annihilating particle 
exceeds the energy range of GLAST, the study of the 
detected sources can be extended to higher energies by 
air Cerenkov telescopes (ACTs) like CANGAROO Q, 
HESS 0, MAGIC H or VERITAS [nj. 

This article is organized as follows. In Sec. HH we give 
an overview over indirect DM searches through gamma 
rays and introduce the gamma-ray spectra that are typ- 
ically expected from the annihilation of DM particles. 
The main characteristics of GLAST, along with a brief 
account of its scientific goals and expected performances, 
are presented in Sec. IIIII In Sec. IIV1 we show the sensi- 
tivity maps for the detection and identification of point 
sources of DM annihilation. We then apply our results 
to mini-spikes around IMBHs in Sec|Vl where we make 
detailed predictions about the detectability of these ob- 
jects. Finally, we present our conclusions in Sec. IVIII 

II. INDIRECT DARK MATTER SEARCHES 
THROUGH GAMMA RAYS 

A theoretically particularly well-motivated type of 
DM candidates are weakly interacting massive particles 
(WIMPs) that appear in various extensions to the stan- 
dard model of particle physics (SM); with masses and 
couplings at the electroweak scale, they would be ther- 
mally produced in the early universe and automatically 
acquire the necessary relic density to account for the DM 
today. Usually, the WIMP appears as the lightest of a 
whole set of new, heavy particles and its decay into SM 
degrees of freedom is protected by a new symmetry. 

The prototype exam ple for a DM candidate of this type 
is the neutralino (see [ll[ for a classic review) that ap- 
pears in most supersymmetric extensions to the SM as 
the lightest supersymmetric particle (LSP) and is given 
by a linear combination of the superpartners of the gauge 
and Higgs fields, 

X = Xi = N U B + N 12 W 3 + N 13 H° + N U H% . (1) 

The neutralino mass is usually several hundred GeV or 
less; for very high Higgsino or Wino fractions, however, 
it can be considerably higher (up to 2.2 TeV in the lat- 
ter case). More recently, considerable attention has also 
turned to models with universal extra dimensions (UED) 
[l2| . where the higher-dimensional degrees of freedom 
manifest themselves as towers of new, heavy states in 
the low-energy, effectively four-dimensional theory. The 
lightest of these Kaluza-Klein particles (LKP) has a mass 
very close to the inverse of the compactification scale R, 
and for ~ 1 TeV it turns out to be aperfect WIMP 
DM candidate just as the neutralino (rj In the 

minimal UED setup [r|, the LKP is given by the 7^, 
the first KK excitation of the photon. For the discussion 



that follows, we will mainly have these DM candidates 
in mind; however, we will also argue that the gamma- 
ray signals we expect in these cases are rather generic, so 
that our conclusions about the prospects for a detection 
of DM annihilation point-sources will remain essentially 
unchanged even for other DM candidates. 

Both LSP and LKP are charged under a Zi symmetry 
(i?-parity and KK parity, respectively), under which all 
SM particles remain neutral. This ensures their stability 
against decay, but at the same time allows the (pair-) 
annihilation into SM particles. The resulting differential 
gamma-ray flux at earth is then given by 



S 7 (£) = 




where the integral over the DM density p is along the line 
of sight for a given angle i\) of observation, (av) is the to- 
tal annihilation rate, m x the dark matter particle's mass 
and E = m x x the photon energy; the sum, finally, runs 
over all possible annihilation channels, with branching 
ratios and photon multiplicities dN^/dx. The left- 
hand side of the above expression for $ 7 (-E) determines 
the spectral form of the signal; it depends only on the un- 
derlying model of particle physics and can be computed 
to a very good accuracy. However, since the dark mat- 
ter distribution p(r) is only poorly known, the same is 
unfortunately not true for the line-of-sight integral that 
appears on the right-hand side of Eq. ([2]). For classical 
targets of indirect DM searches like the galactic center, 
for example, these astrophysical uncertainties result in a 
predicted amplitude for the flux that may vary over sev- 
eral orders of magnitude [l5[ ■ Another problem that indi- 
rect DM searches have to face is that regions of expected 
high DM density usually also contain a sizeable amount 
of ordinary matter which gives rise to more conventional, 
astrophysical sources of gamma-rays. This is particu- 
larly true for the galactic center and even if the observed 
gamma-ray signal from that direction is still not iden- 
tified unambiguously, it is very unlikely that the HESS 
source can mainly be attributed to DM annihilations [16j | . 
A DM component could still be singled out wth GLAST 
in the 10 - 300 GeV energy range, where a DM source can 
provide a good fit of the EGRET data [13, d, but in or- 
der to be conclusive, it is of great importance to provide 
smoking-gun evidence for any claimed indirect detection 
of DM. Such evidence may come from a multi- wavelength 
analysis of the annihilation spectrum (l9| or from pro- 
nounced spectral features that cannot be mimicked by 
astrophysical processes, like line-signals [20( or the very 
sharp cutoffs that are associated with final state radia- 
tion [2l|, [H [H, Hi]. A third possibility is the one that 
we will discuss later, namely the appearance of a great 
number of sources with identical spectra as expected in 
the mini-spike IMBH scenario. 

Let us now investigate in a bit more detail the gamma- 
ray spectra that may result from DM annihilations. In 
principle, there are three different types of contributions 
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FIG. 1: Photon multiplicity for various final states. The 
(blue, magenta, cyan, black,red) solid lines show the (uu/dd, ss, 
cc, bb, tt) quark spectra, while the (green,blue) dotted lines 
give the contributions from (WW, ZZ) gauge bosons. The 
(black) dashed line, finally, is the FSR spectrum from e + e~ 
final states. All spectra are plotted for both m x = 500 GeV 
and m x = f 000 GeV. 




x — E/m x 

FIG. 2: Same as Fig. [T] (for m x = 500 GeV and without 
the e + e~ FSR case), but with spectra rescaled as in ((5J and 
scaling parameters as given in Tab. U The straight dotted 
line represents the extrapolated, asymptotic x~ 15 behaviour 
of the photon spectrum from 66 final states. The thick black 
line shows the (not rescaled) bb spectrum that we will use for 
reference. 



to consider. The dominant source of gamma rays is 
usually associated to secondary photons from quark (or 
gauge boson) fragmentation and decay, mainly through 
the process tt° — > 77. The second contribution arises 
from final state radiation (FSR), where an additional 
photon is emitted from charged particle final states; this 
becomes particularly important for a sizable branching 
ratio into e + e~ pairs [2l|, [22|], though, for m x >1 TeV, 
the same can be true even for W + W~ pairs [2J|. Fi- 
nally, dark matter particles can annihilate directly into 
photons, giving rise to a characteristic line signal at the 
DM particle's mass [2(| ■ However, since the DM has 
to be electrically neutral, these channels (xx ~^ 77 or 
XX ~ * %l) are necessarily loop-suppressed and thus usu- 
ally negligible - at least for current detector resolutions 
and in the absence of efficient enhancement mechanisms 
[25l | . Having in mind the limited reach of GLAST for very 
high energies, we therefore concentrate in the following 
on secondary photons and FSR photons from e + e _ final 
states. 

For the secondary photons, we have adopted the pa- 
rameterizations that were obtained in Ref. [To] ] by run- 
ning the Pythia Monte Carlo code f2f| : 



U-V 7 _ a b+cx+dx 2 +ex 3 

Ijdj O 

ax 



(3) 



where r\ = 2 for the annihilation into WW, ZZ and tt, 
and 1 otherwise. The value of the parameters (a, b, c, d, e) 
for the various annihilation channels / and a set of DM 
masses can be found in Ref. [l5[ . We limit ourselves here 
to show in Figure [1] the photon multiplicities dN^/dx 
per annihilation into quark and gauge boson pairs, for 
two different values of the DM mass, M — 500 GeV and 
M = 1000 GeV. For the bb channel, which we will con- 
sider below, the parameter values are (77, a, b, c, d, e) = 
(1, 1.50, 0.37, -16.05, 18.01, -19.50). 



The e + e~ FSR spectrum, on the other hand, can be 
determined analytically and is, in leading logarithmic or- 
der, given by [22LI23I] 



d/V« _ d((i e+e - 7 )/da; a (x 2 - 2x + 2) 



da; 



In '-^(1-x) 

% 

& 

For the other charged leptons, one has simply to replace 
the electron mass m e in the above expression with the 
corresponding lepton mass, leading to roughly the same 
spectral form as in the e + e~ case. The situation for 
t leptons is actually a bit more complicated since their 
semi-hadronic decay produces secondary photons on top 
of the FSR contribution. The resulting spectrum exhibits 
a power-law behaviour for small photon energies, with a 
spectral index of about —1.3, and a cutoff that is sharper 
than in the case of quarks or gauge bosons, but less pro- 
nounced than for e + e~ pairs [15j. In the following, we 
will not consider this type of spectrum separately. We 
note that even in the LKP case, where gamma-rays from 
t leptons are known to give an unusually important con- 
tribution [12], this does not change our general conclu- 
sions about the determination of the cutoff. 

We are now in a position to make two important ob- 
servations about the gamma-ray spectra that can be ex- 
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uu 


dd 


ss 


cc 


tt 


WW 


ZZ 




1.00 


2.74 


2.74 


1.83 


1.38 


1.36 


3.56 


3.02 


Bf 


1.00 


1.47 


1.47 


1.16 


1.07 


1.05 


1.52 


1.45 



TABLE I: With these scaling parameters, the photon mul- 
tiplicities for the above final states are rescaled according to 
((5}, so as to feature the same asymptotic behaviour and cutoff 
as bb final states (here, we have defined the cutoff as a 10 -3 
drop of the spectrum with respect to the asymptotic a; -1,5 
power law). 
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pected from DM annihilations. First, expressed in the 
variable x = E/m x , they are to a very good approxi- 
mation independent of the DM mass m x . Second, there 
are essentially only two types of spectra. FSR spectra, 
to begin with, show an asymptotic x~~ x power law for 
small x and a characteristic, extremely sharp cutoff at the 
highest energy that is kinematically available. The spec- 
tra of secondary photons, on the other hand, exhibit an 
asymptotic a; -1 ' 5 power law and a much less pronounced 
cutoff. In fact, spectra of this second type are almost 
indistinguishable if the dark matter mass is not known 
by some independent measurement. To make this point 
more clear, let us consider a rescaled version of the pho- 
ton multiplicities for secondary photons, 



dx 



(X) EE Af 



(5) 



Here, Af represents an overall normalization of the am- 
plitude while Bf corresponds to a shift in the cutoff, 
i.e. the DM mass. Since neither of these parameters is 
a priori known, we should use diV^' res /dx, rather than 
dN^ /dx, to compare the different spectra. The scaling 
parameters Af and Bf can then be fixed by demanding 
that all spectra have the same normalization for small 
x and the same cutoff (which we define here, somewhat 
arbitrary, as a 10 -3 drop of the spectrum with respect to 
the asymptotic power law). The resulting rescaled 

spectra are shown in Fig. [2l with scaling parameters as 
given in Table [IJ 

For illustration, we will in the following consider two 
different, rather simplified DM models: in model A, the 
DM mainly annihilates into bb pairs while in model B, we 
fix a branching ratio of 80% into bb and 20% into e + e~ 
final states. The first case corresponds roughly to the 
situation one generally expects for the LSP (at least for 
masses smaller than a few hundred GeV, see, e.g., L 27i\), 
while a large branching ratio into leptons like in model 
B is for example realized in the case of the LKP. From 
the above discussion, however, it should be clear that our 
models A and B are more general in the sense that they 
cover the two basic possibilities for spectral distributions 
that can be expected from DM annihilations. At this 
point, a word of caution is in order. While the charac- 
teristic sharp FSR cutoff can be used to unambiguously 
identify the DM mass to an accuracy that corresponds 
to the energy resolution of the detector, a certain theo- 
retical uncertainty enters into the determination of m x 
when only the rather soft cutoff from secondary photons 
is observed. If, e.g., DM mainly annihilates into light 
quarks or gauge bosons, one would not be able to tell the 
difference from our case A, but underestimate the actual 
DM mass by up to 50% (see Table U). 

Let us conclude this Section by stressing once more 
that the detection of many sources with the same spec- 
trum, and in particular with the same cutoff, would pro- 
vide a smoking-gun evidence for DM annihilations - even 
if the cutoff is not very pronounced, like in model A, and 
would therefore not provide conclusive evidence when ob- 



Peak Effective Area (1-10 GeV) 


>8000 cm 2 


Energy Resolution 10 GeV on-axis 


<10% 


Energy Resolution 10-300 GeV on-axis 


<20% 


Energy Resolution 10-300 GeV >60° 


<6% 


PSF 68% 10 GeV on-axis 


< 0.15° 


Field of view 


2.4 sr 


Point source sensitivity (>100MeV,5y) 


< 6 x 10" 9 cm- 2 s _1 


Source Location 


< 0.5 arcmin 



TABLE II: GLAST LAT minimum performances, adapted 
from Ref. El 



served from a single source alone. 



III. GLAST 

The Gamma-ray Large Area Space Telescope 
(GLAST) is a next generation gamma ray observatory 
due for launch in Fall 2007. The main scientific objects 
are the study of all gamma ray sources such as blazars, 
gamma-ray bursts, supernova remnants, pulsars, diffuse 
radiation, and unidentified high-energy sources. Indeed, 
GLAST might be called the "Hubble Telescope" of 
gamma-ray astronomy as it will be able to observe 
AGN sources up to z ~ 4 and beyond, if such objects 
actually existed at such early times in the universe. 
Extrapolation from EGRET AGN detections shows that 
about 5,000 AGN sources will be detected in a 2 years 
cumulative scanning mode observation by GLAST, as 
compared to the 85 that have been observed by EGRET 
in a similar time interval. The GLAST collaboration 
has also a strong interest in DM detection and several 
works in the literature investigate possible scenarios 
and the corresponding detection capabilities (see e.g. 
Ref. [28| and the aforementioned reviews [H, El for 
further details). 

The primary instrument on board is the Large Area 
Telescope (LAT), a pair- production telescope designed 
to allow observation of celestial sources in the range 
20 MeV-300 GeV (2^; LAT was developed by an interna- 
tional collaboration of particle physics and astrophysics 
communities from 26 institutions in the United States, 
Italy, Japan, France and Germany. The LAT is an 4x4 
array of identical towers, each formed by 

• Si-strip Tracker Detectors and converters arranged 
in 18 XY tracking planes for the measurement of 
the photon direction, 

• Segmented array of CsI(Tl) crystals for the mea- 
surement the photon energy, 

• Segmented Anticoincidence Detector (ACD). 

During normal operation the LAT will continuously 
scan the whole sky, obtaining a complete coverage about 
every 3 hours due to the great field of view (~ 2 sr). The 
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uniform sky coverage, in addition to the large effective 
area and the good angular resolution should ensure many 
advances in the study of all celestial sources of interest. 
In Tab.HIl we summarize from [301 ] the LAT performances 
more relevant in our context. 



IV. DARK MATTER SENSITIVITY MAPS 

In order to study the LAT sensitivity for DM annihila- 
tion signals, wc perform a simulation of the gamma-ray 
sky as it will be observed by this instrument, based on 
a parametrization of the instrument response. This al- 
lows us to rapidly simulate several months of data-taking 
in a relatively short time, without having to deal with a 
complex MonteCarlo simulation of the physical processes 
involved. For additional information about GLAST LAT 
simulations, see, for example, Ref. (3l| . 

All simulations performed for this study correspond to 
2 months of scanning (i.e. not pointed) LAT observations. 
The time duration was selected such as to allow for a rel- 
atively fast discovery of a DM signal; the trade-off is a 
sizeable decrease in the number of resolved and correctly 



identified sources. No other (non-DM) point sources were 
simulated and source spatial coincidence was not consid- 
ered, in order to simplify the analysis and to allow for an 
automated analysis process. 

The extragalactic diffuse background was simulated by 
extending the EGRET observations to the LAT energy 
range [321 ] . The galactic diffuse background model is im- 
plemented in the framework of the GALPROP code for 
cosmic-ray propagation and incorporates up-to-date sur- 
veys of the interstellar medium, as well as current mod- 
els for the interstellar radiation field, updated production 
functions and inverse scattering calculations [331 ] . As mo- 
tivated in Sec. II, we will focus on DM candidates mainly 
annihilating into bb, and discuss below how the detection 
of deviations from the corresponding gamma-ray spec- 
trum may provide useful insights on the nature of the 
DM particles. 

In the first part of our analysis, we investigated the 
intensity required to resolve a point-source of DM an- 
nihilations at a 5cr significance level. We divided the 
sky into regions of about 10 degrees in radius, and in 
each region we placed one DM source. Then, we consid- 
ered each source separately and let the flux intensity vary 




FIG. 4: GLAST sensitivity map for the detection of point sources of Dark Matter annihilations, i.e. full-sky map in galactic 
coordinates of the minimum flux above 20 MeV, in units of [ph m _2 s _1 ], that is required for a 5<r detection of an annihilation 
spectrum, assuming a DM particle with mass m x — 150 GeV annihilating into bb (note, however, that the map does not depend 
very sensitively on DM properties). The map is relative to a 2 months observation period; for longer observation times, fluxes 

— 1/2 

scale approximately as t ob J . For reference, we also show the simulated gamma ray sky. 




FIG. 5: GLAST sensitivity map for the identification of point sources of Dark Matter annihilation, i.e. full-sky map in galactic 
coordinates of the minimum flux above 20 MeV, in units of [ph m _2 s _1 ], necessary to discriminate the annihilation spectrum 
from an ordinary astrophysical source, assuming a DM particle with mass m x — 150 GeV annihilating into bb, after a 2 months 
observation period (see text for a further discussion). For reference, we also show the simulated gamma ray sky. 



from 10 -4 to 2 x 10 -3 ph m _2 s _1 . For each intensity, we 
calculated the significance of the observed signal, given 
the local background counts, with a maximum likelihood 
analysis assuming Poisson statistics. By estimating the 
minimum flux required to discriminate the DM source 



from the background at a 5cr level on a grid of points 
uniformly distributed over the sky, we have obtained the 
sensitivity map shown in the top panel of Fig. [4] (where 
we adopted a DM particle mass m x = 150 GeV). The 
sensitivity appears to depend significantly on the Galac- 
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tic longitude only along the Galactic disk, as expected. 
At high galactic latitudes a source as faint as 2 x 10~ 4 
ph m _2 s~ 1 is resolved, while close to the galactic center 
a minimum flux of 1.5 x ph m _2 s _1 is required. 

In a similar way, we can also produce a sensitivity 
map for a reliable identification of the DM spectral cut- 
off. To do so, we compare for each observed spectrum 
the likelihood of a plain power law, with that of a DM- 
like spectrum, and check whether the latter hypothesis is 
favoured. The result of this analysis is shown in Fig. 
Inside the galactic plane (\b\ < 5) the minimum intensity 
required to pinpoint the high energy cutoff is far greater 
than the intensity required for a 5a detection. This im- 
plies that there will be some sources that can correctly 
be identified as a DM signal only after observation times 
longer than the 2 months considered here. 

All figures quoted in this work were calculated by 
taking into account in the analysis only detected pho- 
tons with energies between 100 MeV and 200 GeV. We 
limit ourselves to this energy range in order to avoid the 
(present) uncertainties in the LAT response functions at 
very low (down to ~ 20 MeV) and high energies (up to 
500 GeV). On the other hand, we do take into account all 
major effects such as orbiting, rocking and dead times. 
While our sensitivity maps are in principle still subject to 
further revisions as the knowledge of the LAT response 
improves, we are therefore confident that they will re- 
ceive only minor corrections (improving the detection ef- 
ficiency) . 

The time scale of 2 months adopted here corresponds 
roughly to the periodicity of the LAT orbit given an in- 
clination of ~ 28.5° (~ 55 days). Assuming that the 
detection statistics are dominated by diffuse counts and 
not by stochastic fluctuation, it is possible to build de- 
tection maps for longer observation times t b s by sim- 

— 1/2 

ply scaling down the minimum flux with t ob J . For in- 
stance, by increasing the observation time by a factor 4, 
t bs = 8 months, the minimum fluxes shown in the maps 
are roughly halved. For the identification maps, on the 
other hand, the situation is more complicated. It is in 
fact not trivial to scale minimum fluxes for the determi- 
nation of the DM mass to longer observation time, since 
the analysis is based on a comparison of likelihood be- 
tween models, in a regime where both background counts 
and low statistics play an important role. Although a 
scaling t^ 2 appears as a conservative estimate in this 
case, we feel that this issue deserves a dedicated study. 

The sensitivity map for the detection of DM point 
sources is rather insensitive to the DM particle mass, 
since the flux is dominated by the low energy tail of the 
annihilation spectrum. However, things are different for 
the identification map, at least above ~ 200GeV. In fact, 
the determination of a high energy cut-off strongly de- 
pends on the possibility of obtaining good quality data 
(in terms of angular and energy resolution), for photons 
above this energy. A detailed study of high-mass mod- 
els will thus be possible only with a reliable simulation 
of the instrument response at energies above the range 



considered, although it should be possible to identify its 
DM origin, through the rather distinctive slope of the 
annihilation spectrum (see discussion below). Finally, 
we note that neither the detection nor the identification 
of DM point sources depends significantly on the annihi- 
lation channels; this situation would, in principle, change 
for longer observation times and DM masses below a few 
hundred GeV, when the pronounced sharp cutoff for FSR 
photons from e + e~ final states becomes clearly visible, 
thereby enhancing the prospects for the identification of 
model B type DM. 

We stress that the sensitivity maps in Fig. [^] and [3| 
represent a powerful tool to assess the observability of any 
population of point sources of DM annihilations, without 
having to perform a case-by-case study (simulation plus 
analysis) of individual sources. 

V. DETECTABILITY OF MINI-SPIKES 

We now apply the results obtained in the first part 
to a specific astrophysical scenario, namely to DM mini- 
spikes around Intermediate Mass Black Holes (IMBHs), 
and study the prospects for detecting these objects with 
GLAST. In this section, following closely Ref. (34|, we 
will first review the basics of the mini-spike scenario, and 
derive, for a Milky- Way like galaxy, their number, pro- 
file and total luminosity in terms of DM annihilations. 
We will then combine this information with the sensitiv- 
ity maps shown above to determine the prospects for a 
detection with GLAST. 



A. The Mini-spikes scenario 

The effect of the formation of a central object on the 
surrounding distribution of matter has been investigated 
in Refs. [H, [H, H3, HI] and for the first time in the 
framework of DM annihilations in Ref. [3^| . It was shown 
that the adiabatic growth of a massive object at the center 
of a power-law distribution of DM with index 7, induces 
a redistribution of matter into a new power-law (dubbed 
"spike") with index 



This formula is valid over a region of size r sp w 0.2r^, 
where tbh is the radius of gravitational influence of the 
black hole (see below for further details). The adiabatic- 
ity of the black hole's growth is in particular valid for the 
SMBH at the Galactic center. A critical assessment of the 
formation and survival of the central spike, over cosmo- 
logical timescales, is presented in Refs. EH (see also 
references therein). We limit ourselves here to note that 
adiabatic spikes are rather fragile structures, that require 
fine-tuned conditions to form at the center of galactic ha- 
los [13], and that can be easily destroyed by dynamical 
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processes such as major mergers [43[ and gravitational 
scattering off stars jlolEj . 

It was recently shown that a p oc r~ 3 / 2 DM overden- 
sity can be predicted in any halo at the center of any 
galaxy old enough to have grown a power-law de nsity 
cusp in the stars via the Bahcall-Wolf mechanism [4a |. 
Collisional generation of these DM "crests" (Collision- 
ally REgenerated STtructures) was demonstrated even 
in the extreme case where the DM density was lowered 
by slingshot ejection from a binary super-massive black 
hole. However, the enhancement of the annihilation sig- 
nal from a DM crest is typically much smaller than for 
adiabatic spikes (45l |. 

Although it is unlikely that a spike may survive around 
the Super-massive Black Hole at the Galactic center, 
they can evolve unperturbed around Intermediate Mass 
Black Holes (IMBHs), i.e. wandering BHs with mass 
10 2 < M/Mq < 10 6 . Scenarios that seek to explain 
the properties of the observed super-massive black holes 
population result, in fact, in the prediction of a large pop- 
ulation of IMBHs. The number and properties of mini- 
spikes around IMBHs have been discussed in Ref. [34| . 
where two different formation scenarios have been inves- 
tigated. In the first one, IMBHs form in rare, overdense 
regions at high redshift, z ~ 20, as remnants of Popula- 
tion III stars, and have a characteristic mass-scale of a 
few 1O 2 M [46| (the formation of mild mini- spikes around 
these objects was also discussed in Ref. [13, IH, 1491). 
In this scenario, these black holes serve as the seeds for 
the growth of super-massive black holes found in galactic 
spheroids [5(|. We will not further consider this scenario, 
because of the many uncertainties associated to the ini- 
tial and final mass function of IMBHs, and the lack of a 
solid prescription for the initial DM distribution and BH 
mass grow due to gas accretion as these objects move in 
galactic halos. 

We will instead focus on the second scenario discussed 
in Ref. [34| , based on the proposal of Ref. [5l| , which is 
representative of a class of models in which black holes 
originate from massive objects formed directly during 
the collapse of primordial gas in early- forming halos. In 
practice, during the virialization and collapse of the first 
halos, the gas cools and forms pressure-supported disks 
at the centers of halos that are sufficiently massive for 
molecular hydrogen cooling to be efficient. Gravitational 
instabilities in the disk lead to an effective viscosity that 
transfers mass inward and angular momentum outward 
until the first generation of stars heats the disk and ter- 
minates this process. In this case, the characteristic mass 
of the black hole forming in a halo of virial mass M v is 
given by [5l[ 

M hh =3 .8xlO 4 M (^)( I 4) 3/2 

\^10 7 A% J \ 18 / ^lOMyry ' \' J 

where / is the fraction of the total baryonic mass in the 
halo that has fallen into the disk, z is the redshift of 



formation, k is the baryon fraction that lost its angular 
momentum on a time t ev , and t ev is the timcscalc for the 
evolution of the first generation of stars [5l[ . The distri- 
bution of black hole masses is a log-normal distribution 
with a mean value given by the characteristic mass above 
and a standard deviation OM bh = 0.9. 

The "initial" DM host mini-halo (that is, the DM dis- 
tribution prior to black hole formation) can be well ap- 
proximated by a Navarro, Frenk, and White (NFW) pro- 
file m 

*) _„,(£)-'(, + £)-'. (8) 

The normalization constant po, and the scale radius r s , 
can be expressed in terms of the virial mass of the halo 
at the time when the IMBH formed, M V \ T , and the virial 
concentration parameter, c V i r . See Ref. [3J] for further 
details. Alternatively, we could have chosen the more re- 
cent parametrization proposed by Navarro et al. (53j (see 
also Refs. [54|, |55[). However, this profile implies mod- 
ifications at scales smaller than those we are interested 
in, where the profile is anyway modified by the presence 
of the IMBH, hence leading to a negligible change in the 
predicted fluxes. 

We assume that the black holes form over a timescale 
long enough to guarantee adiabaticity, but short com- 
pared to the cosmological evolution of the host halo (see 
Ref. [HI for further details). In fact, the condition of 
"adiabaticity" , fundamental to grow "mini-spikes" , re- 
quires that the formation time of the black hole is much 
larger than the dynamical timescale at a distance from 
the black hole, where is the radius of the sphere of 
gravitational influence of the black hole, ~ GM^/a 2 , 
and a is the velocity dispersion of DM particles at r^. In 
practice, is estimated by solving the implicit equation 

M{< r h ) = f " p(r)r 2 dr = 2 M bh . (9) 
Jo 

Inside a region r sp ss 0.2r^ [56], and assuming a NFW 
profile as in Eq. ([8]). the spike profile is given by 

Psp = p(r S p)[ — ) . (10) 

\ r sp / 

Such a steep power-law leads to a density profile that di- 
verges at small radii. Although the Schwarzschild radius 
of the black hole provides already a physical lower cut- 
off, the profile already 'saturates' at larger radii, where 
the DM annihilation timescale becomes shorter than the 
time t sp elapsed since the formation of the spike. This 
radius ru m can be expressed implicitly as 

Ps P (r cu t) = m x /av (t - t sp ) , (11) 

where av is the annihilation cross section times the rel- 
ative velocity of the DM particle. When performing our 
integrations, we thus stop at a "cut" radius defined as 

r cut (to, av) = Max [4R s ,r lim (to, av)] , (12) 
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resolved sources 




FIG. 6: Number of realizations Nr where a number Ns of mini-spikes is observed (empty histogram) and identified as DM 
sources (red dashed), assuming m x = 150 GeV. 





Accepted realizations 


Resolved mini-spikes 


Identified mini-spikes 


Standard scenario 


159/200 


22 ±8 


19 ±7 


Scaled £ = 0.2 


198/200 


4±2 


4±2 


Scaled £ = 5 


13/200 


52 ± 14 


45 ± 12 



TABLE III: Statistics of the number accepted realizations (i.e. those not violating the EGRET constrain), resolved sources 
(i.e. those detected by GLAST within 2 months at 5cr) and identified sources (i.e. those that can be discriminated by ordinary 
astrophysical sources), for M=150 GeV, and 2 different values of the rescaling factor £ that accounts for a variation of the 
astrophysical and particle physics parameters in the mini-spike model. 



where R s = 2.95 km M^/Mq is the Schwarzschild radius 
of the IMBH. 

The flux of gamma-rays from a mini-spike around an 
IMBH can be expressed as 34] 



dN 



AE VlCr 26 cm 3 /s 
kpc 



)( 



in. 



100GeV/ 

2 



pc 



10 2 GeVcm" 3 
r cut (m, eru) 



10" 3 pc 



(13) 



where $ = 9 x 10 _10 cm~ 2 s _1 and D is the IMBH dis- 
tance to the Earth. Note that the annihilation flux in 
this case does not depend, as one may naively expect, on 
the ratio (<ti>/to 2 ), because r cut is itself a function of the 
particle mass and cross section. It is easy to show that in 
this case the (integrated) gamma-ray flux is proportional 
to ~ (av) 2 / 7 m~ 9/J ( given that dN/dx is almost indepen- 
dent of m x ), which means that - due to the "saturation" 
of the DM profile described above - the predicted fluxes 



are much less sensitive to the particle physics parameters 
than in the usual case. 

In order to assess the detectability of mini-spikes, we 
have followed the procedure outlined in Refs. [H, H3] 
to build mock catalogs of sources. As detailed in the 
next Section, we have then performed our analysis on 200 
statistical realizations of a Milky- Way like halo, studied 
the detectability in each realization, and then performed 
a statistical analysis of the results. 



B. Prospects for detection with GLAST 

We now turn to the prospects for detecting mini-spikes 
with GLAST. To this aim, we simply have to compare 
our mock catalogs of sources, described in more detail in 
the previous section, with the sensitivity maps of Figs.0] 
and 03 If the flux from a given source is greater than 
the minimal flux required to grant a 5cr detection we will 
say that the source is resolved; if it is also greater than 
the threshold for the discrimination from astrophysical 
sources, we will say the source is identified. Let us point 
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FIG. 7: Number of realizations Nu where a number Ns of mini-spikes is observed (empty histogram) and identified as DM 
sources (red dashed), assuming m x = 150 GeV, and a rescaling factor £ = 0.2 (left) and £ = 5 (right). 



out here that, obviously, the discussion on the detectabil- 
ity of mini-spikes can only be performed in a statistical 
sense, since we do not know in what realization of the 
simulated Milky Way halo we live in. 

The result of the above procedure is shown in Fig. [6l 
where we show the number Nr of realizations with Ns 
resolved (and identified) mini-spikes. Let us stress that 
some of the realizations actually contain sources that 
would shine brighter than the brightest EGRET uniden- 
tified source outside the galactic plane (3EG J1835p5918, 
[58|). In order not to violate the observational con- 
straints, we therefore reject any realization where the flux 
of at least one source exceeds $ max = 2 x 10 -2 ph m~ 2 
s" 1 above 20 MeV; these rejected realizations are not 
included in Fig. [6] and will not be considered in the fol- 
lowing, either. Incidentally, it is interesting to note that 
the spectral index reported for the brightest unidenti- 
fied EGRET source is n = 1.69 ± 0.07, quite close to 
what is expected for the low energy tail of the annihila- 
tion spectrum (EGRET was sensitive to energies up to 
30 GeV). At lower energies, COMPTEL (in the range 
0.75-30 MeV) did not resolve this source from the steep 
blazar QS01739+522 [H. 

We can see from Fig. [5] that, in most realizations and 
adopting a particle mass m x = 150 GeV, we can expect 
a large number of sources to be detectable; Tab. IIII1 fur- 
thermore, shows that only a small fraction of all the re- 
alizations is in violation of the observational constraints 
set by EGRET. 

The expected flux from any given mini-spike does, 
however, depend on various particle physics and as- 
trophysical parameters. For what concerns the parti- 
cle physics parameters, we have already stressed that 
&i(E) cx (ctu) 2 / 7 to x 9 ^ 7 . As for the astrophysical uncer- 
tainties, we have shown in Eq. [7] the dependence of the 
BH mass on the formation model parameters. Without 
entering here into a detailed discussion of the parameter 
space of the model, we incorporate the effect of varying 
all parameters, by introducing a global rescaling param- 



eter £, building new realizations with fluxes 

=£*«(£?,£>) , (14) 

where $i are the fluxes obtained in the reference stan- 
dard scenario discussed above, which then corresponds 
to the case £ = 1. For high values of £, we expect a 
lower number of acceptable sources, since many of them 
will violate the EGRET constraint, while lowering £ will 
make a larger number of realizations compatible with the 
same bound. We show the results for the cases £ = 0.2 
and 5, respectively, in Fig. and summarize the number 
of accepted realizations, resolved sources and identified 
sources, in Tab. IIIII 

VI. TOWARDS THE NATURE OF DARK 
MATTER 

In this Section, we analyze in more detail some aspects 
of the detected sources. As a first step, it is important to 
verify under which conditions the LAT is able to detect 
the cutoff for a DM candidate like, e.g., in our model A, 
i.e. under which conditions the fit obtained by adopting 
the parametrization of Eq. ([3]) (with only m x and the 
absolute normalization as free parameters) has a greater 
likelihood than an ordinary power-law fit (where the free 
parameters are the overall normalization and the spectral 
slope). We find that the highest DM mass m x that still 
gives a spectrum recognizable as DM-induced is close to 
1 TeV for the brightest sources. The efficiency in the m x 
estimate, however, is strongly dependent on the mass it- 
self. For a simulated source with m x = 150 GeV, e.g., we 
can estimate the mass of the DM particle with a precision 
better than 25 GeV. For higher masses the precision on 
the determination of the DM mass depends strongly on 
the possibility to extend the instrument's energy range 
to photon energies above 200 GeV; an educated guess on 
the instrument's response suggests that a mass of 1 TeV 
can be reconstructed with a precision of 10% for the 
brightest sources. A detailed analysis would warrant a 
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FIG. 8: Examples of spectral fits of simulated DM point sources of intensity <j>, for different values of m x and different 
annihilation channels. Solid lines are fits obtained under the assumption of annihilation to bb. For each model we also give the 
significance of the detection. Points with error bars are photon counts from the simulated observation. 



dedicated study, also because the exact role of system- 
atics and stochastical fluctuation is at present not clear. 
Mass spectroscopy, in any case, should be postponed un- 
til a better knowledge of the instrument's response at 
high energies and data from longer observation periods 
become available. 

Turning to some illustrative examples of simulated 
sources, we show in Fig. [5] the situation of a model A - 
type source, at a position (l,b)=(0,25) in galactic coordi- 
nates, with a mass m x = 150 GeV and a flux <j> = 2 x 10 -3 
ph m~ 2 s _1 . At this location, we expect a moderate dif- 
fuse background contribution; it was chosen such as to 
demonstrate that a source corresponding to EGRET's 
faintest detected source is easily identified. In the same 
figure, we also present a source with m x = 1 TeV and 
(f> = 2 x 10~ 2 ph m~ 2 s^" 1 . Here, we have picked a position 
in the galactic plane, at (l,b)=(50,0), to show that even 
in presence of an intense diffuse gamma-ray background, 
a cutoff at 1 TeV is well within the reach of our analysis 
- at least for the brightest sources in our realizations. 
In the right panel of Fig. [5J finally, we show an example 
spectrum for a bright (<j> = 2 x 10~ 2 ph m~ 2 s _1 ) model 
B - type DM source with m x = 150 GeV. The source 
was simulated at galactic coordinates (l,b)=(0,50), where 
background counts are negligible; the superimposed fit 
was computed with the usual Eq. ([3]), i.e. without taking 
into account the leptonic contribution. At high energies, 
one clearly notices a photon excess due to the leptonic 
component; however, the photon counts are nonetheless 
rather small and therefore fluctuations play a major role. 
As observation time increases, a leptonic component as 
in model B would be recognizable as such even for fainter 
sources; for the brighter sources, on the other hand, we 
expect the characteristic, sharply pronounced cutoff to 
become ever more visible. For a more detailed discus- 
sion, a dedicated analysis would be necessary; we leave 
this interesting issue open for future studies. 



To conclude this Section, let us briefly return to the 
question of how to discriminate a DM from a mere as- 
trophysical origin of gamma-ray point sources, which is 
certainly mostly based on the detection of the high en- 
ergy cutoff at the mass of the DM particle. The GLAST 
LAT is expected to observe more than a hundred of Ac- 
tive Galactic Nuclei (AGNs) in the period we selected 
for this analysis (2 months), some of which would, due to 
EBL attenuation [60j , feature a spectral signature similar 
to that of DM annihilations. However, a population of 
DM sources would exhibit distinctive features that should 
allow to discrminate the from AGNs, such as the low en- 
ergy spectral index of ~ 1.5 and the absence of variability. 
Furthermore, a population of DM sources should appear 
as a distinctive excess in the Fazio-Stecker plot (see for 
example [6(|) that can be obtained from the LAT AGN 
catalog. 



VII. DISCUSSION AND CONCLUSIONS 

The upcoming space telescope GLAST will perform 
full-sky observations of the gamma-ray sky, at energies 
that are of particular interest for indirect DM studies. In 
this paper, we have performed a 2 months simulation of 
the gamma-ray sky as observed by the LAT instrument 
on board the satellite, and, assuming a gamma-ray spec- 
trum typical for DM annihilations, we have estimated the 
minimum flux, at all galactic latitudes and longitudes, re- 
quired to detect the source at 5cr above the background, 
and to identify it, i.e. to discriminate it against ordinary 
astrophysical sources through the analysis of its spectral 
features. 

We have used this information to produce the sensi- 
tivity maps in Figs. [5] and which represent the main 
result of this paper. The first one is a sensitivity map 
for the detection of point sources of DM annihilation, the 



12 



second a sensitivity map for the identification of its DM 
origin, and determination of the mass thanks to a fit of 
the high-energy cut-off. It should actually be possible 
to identify the DM origin of a population of gamma-ray 
sources also in absence of a precise determination of the 
high-energy cut-off, thanks to the distinctive value of the 
spectral slope in the low-energy tail, as discussed in Sec- 
tion II, the presence of an excess in the Fazio-Stecker 
plot, and absence of variability, as argued in Section VII. 
However, the precise determination of the high-energy 
cut-off would provide the opportunity of making an even 
stronger case for the 'exotic' origin of the gamma-ray 
signal, since astrophysical sources can hardly mimic DM 
spectral features (see also the recent paper by Baltz et 
al. [fUl). 

These maps represent a powerful tool to rapidly, and 
efficiently, investigate the detectability of any population 
of point sources of DM annihilation, such as small scale 
clumps [62|, l63l 164. l6a Irjq. 16711 . compact DM structures 
like Spikes [pTpJill I4TI0I or Crests [H[, and DM 
mini-spikes [34 |47| All these objects would appear 
as pointlike sources with GLAST, if their angular size is 
smaller than the angular resolution of the instrument, i.e. 
~ 0.1°. Objects of physical size R p would thus appear 
pointlike when at distance D > 5.7 x 10 2 i? p , a constraint 
often satisfied by the objects listed above. 

In order to prove the effectiveness of the method, we 
have applied our results to the mini-spikes scenario dis- 
cussed in Ref. [34], consisting of a population of ~ 100 



DM overdensities, dubbed mini-spikes, around Interme- 
diate Mass Black Holes. We were thereby able to show 
that a large number of these objects can be detected and 
identified with GLAST, if they exist, while null searches 
would place extremely stringent constraints on the whole 
scenario. Although a set of fiducial astrophysical param- 
eters was chosen for the mini-spike scenario, in order to 
perform this study, the analysis can easily, and rapidly, 
be generalized to any other set of parameters. 

Finally, we have shown that, in case of detection, 
GLAST can provide useful information on the nature of 
the annihilating DM particles: first, with a relatively ac- 
curate determination of their mass, and second, with the 
identification of some specific spectral features, like a pos- 
sible hard spectrum contribution near the cut-off at the 
DM mass, due to the annihilation to charged leptons. 
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